Spin-orbital gap of multiorbital antiferromagnet 



o 
o 

(N 



Hiroaki Onishi and Takashi Hotta 
Advanced Science Research Center, Japan Atomic Energy Agency, Tokai, Ibaraki 319-1195, Japan 

(Dated: July 24, 2006) 

In order to discuss the spin-gap formation in a multiorbital system, we analyze an eg-orbital 
Hubbard model on a geometrically frustrated zigzag chain by using a density-matrix renormalization 
group method. Due to the appearance of a ferro-orbital arrangement, the system is regarded as a 
one-orbital system, while the degree of spin frustration is controlled by the spatial anisotropy of the 
orbital. In the region of strong spin frustration, we observe a finite energy gap between ground and 
first-excited states, which should be called a spin-orbital gap. The physical meaning is clarified by an 
effective Heisenberg spin model including correctly the effect of the orbital arrangement influenced 
by the spin excitation. 
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PACS numbers: 75.45.-l-j, 71.10.Fd, 75.30.Et, 75.40.Mg 



I. INTRODUCTION 

Exotic properties of low-dimensional antiferromagnets 
have continued to attract much attention in the re- 
search field of condensed matter physics, since strong 
quantum spin fluctuation destroys classical Neel order- 
ing even at zero temperature, while it leads to vari- 
ous types of quantum disordered statesj^ In fact, sub- 
stantial advances in material synthesis and experimen- 
tal probes have made it possible to gain deep insights 
into strange spin-liquid phases and quantum phase tran- 
sitions in realistic low-dimensional compounds. In par- 
ticular, a class of spin-gapped antiferromagnets has pro- 
vided rich macroscopic quantum phenomena such as 
spin-Peierls transition and impurity-induced magnetic 
ordering in CuGeOajSiMi and Bose-Einstein condensa- 
tion of magnons in TlCuCla iS*2i& 

Microscopic aspects of the spin-gapped phase in low- 
dimensional antiferromagnets have been extensively dis- 
cussed on the basis of the antiferromagnetic (AFM) 
Heisenberg model. In general, a spin gap is defined as 
the energy difference between singlet ground and triplet 
first-excited states. The existence of the finite spin gap 
can be understood by the concept of a valence-bond 
state. One typical example is a spin-1/2 Heisenberg 
zigzag chain ifiiifiiiiiiSiii In fact, the combined effects of 
quantum fiuctuation and geometrical frustration yield a 
transition from a gapless spin-liquid phase to a gapped 
dimer phase. In the gapped dimer phase, a valence bond 
of neighboring spins is stabilized, but the correlations 
among the valence bonds are significantly suppressed due 
to the effect of spin frustration, implying that finite en- 
ergy is required to excite the valence bond from singlet 
to triplet. 

In experiments, as for a model material of the zigzag 
spin chain, several compounds have been reported, such 
as SrCu02fi^*i^ Cu(ampy)Br2,i6 and (N2H5)CuCl3iLia 
However, a candidate material with a significant spin gap 
has not been found, since the frustration effect is much 
weak in these materials. In such a circumstance, in order 
to realize the spin-gapped dimer phase in a region of the 
strong spin frustration, it is important to explore the 



way to control the degree of the spin frustration. In this 
context, it is expected that orbital degree of freedom is a 
key ingredient to cause highly non-uniform spin exchange 
interactions. Indeed, since d- and /-electron orbitals in 
recently focused materials are spatially anisotropic, the 
sign and magnitude of spin exchange interactions depend 
on the relative positions of ions and occupied orbitals. In 
such a case, we envision that orbital ordering occurs so as 
to minimize the magnetic energy of the spin system. In 
fact, we have found that orbital ordering actually relaxes 
the frustration effect in the exchange interaction due to 
the spatial anisotropy of relevant orbital, even though 
the lattice originally exhibits a geometrically frustrated 
structure. ^^■^''■^^ 

After the occurrence of orbital ordering, one may 
naively expect that the system is described by an ef- 
fective spin system on the frozen background of orbital 
arrangement. It is true that we can intuitively under- 
stand what kind of magnetic correlation develops in the 
static orbital ordering. However, since spin and orbital 
degrees of freedom strongly interact with each other at 
the same energy scale, spin dynamics is influenced by or- 
bital one in principle, unless orbital degree of freedom is 
completely quenched. Such interplay of spin and orbital 
dynamics has been discussed for spin-orbital modelsi^SiS^ 
Here, we point out that even though we observe a spin 
gap, this would not be a pure spin excitation but in- 
cludes the changes of both spin and orbital states. The 
spin gap is definitely given by the energy difference be- 
tween singlet ground and triplet first-excited states, but 
the orbital arrangement in the ground and excited states 
can be different. We believe that it is an intriguing issue 
to clarify how the spin excitation is described in such a 
variable background of orbital arrangement. 

In this paper, we investigate the spin gap of an Cg- 
orbital Hubbard model on a zigzag chain by numerical 
techniques. Due to the appearance of a ferro-orbital (FO) 
arrangement, the system can be always regarded as a one- 
orbital system irrespective of the level splitting. When 
two orbitals are degenerate, the zigzag chain is decoupled 
to a double-chain spin system to suppress the spin frus- 
tration due to the orbital anisotropy. On the other hand. 
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FIG. 1: (a) Lattice arrangement and site numbering of A'^-site 
zigzag chain, (b) Schematic view of level splitting between Cg 
orbitals. 



the orbital anisotropy is controlled by the level splitting, 
leading to the revival of the spin frustration and a finite 
spin-excitation gap. However, the orbital arrangement is 
found to change according to the spin excitation, suggest- 
ing an important concept of dynamical interplay between 
spin and orbital degrees of freedom. 

The organization of this paper is as follows. In Sec. II, 
we describe the model Hamiltonian and numerical tech- 
niques. In Sec. Ill, we show the results of the spin and 
orbital structures in the ground state. Taking account 
of the orbital arrangement, we will introduce an effective 
Heisenberg model. Then, in Sec. IV, we discuss the spin- 
excited states based on the effective model. Finally, the 
paper is summarized in Sec. V. 



II. MODEL AND METHOD 

Let us consider eg orbitals on a zigzag chain with N 
sites. The lattice is composed of equilateral triangles 
and located in the xy plane [see Fig. 1(a)]. Note that 
the zigzag chain is considered as a two-chain system con- 
nected by a zigzag path. We consider the case of one 
electron per site (quarter filling), which is correspond- 
ing to, for example, low-spin state of Ni'^^ ion with 3(f 
configuration. 

When the Hund's rule interaction is small, the ground 
state is found to be paramagnetiCfi2i2S which is relevant 
to a geometrically frustrated antiferromagnet. Here, we 
study the property of the paramagnetic phase and simply 
ignore the Hund's rule interaction, but the results do not 
change qualitatively for quarter filling case. Then, the 
Hamiltonian is given by 

H ^ Y. Cr'dlA+'^r'. - (A/2) Y^{p,a - Pib) 

i,a,r,r',cr i 

+ UYpir^Piri+U'YpiaPib, (1) 



where diao- (diba-) is the annihilation operator for an elec- 
tron with spin a in the 3z^—r^ [x^—y^) orbital at site i, 

PiTa^d\^^dira, Pir^J^a Pircr , a is the vcctor Connecting 
adjacent sites, and i^^, is the hopping amplitude between 
T and r' orbitals along the a direction. 

The hopping amplitudes are evaluated from the over- 
lap integral between gg orbitals in adjacent sitesr^i2SiS& 
which are given by 
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FIG. 2: (a) Original zigzag chain including two orbitals at 
each site, (b) Superblock configuration, in which each orbital 
is regarded as a site. The dotted line denotes the path for the 
finite-system algorithm. 



for the X direction, 

C = V4, C - tZ = V3i/8, tj^b = 3V16 (3) 

for the u direction, and i"~,^=i"^/. Thus, the hopping 
amplitudes depend on the bond direction and occupied 
orbitals. Hereafter, t is taken as the energy unit. 

The second term of H denotes the crystalline elec- 
tric field (CEF) potentials. Although the CEF potential 
should be evaluated by the sum of electrostatic potentials 
from ligand ions surrounding the transition metal ionjS^ 
here we consider the level splitting A between 3z^— r^ 
and x'^—y^ orbitals due to the tetragonal CEF effect, as 
shown in Fig. 1(b). In the third and fourth terms, U and 
U' indicate intraorbital and interorbital Coulomb interac- 
tions, respectively. For the local interaction parameters, 
the relation U=U' holds due to the rotational invariance 
in the local orbital space, when we ignore the Hund's rule 
coupling and pair hopping parameters .—■ 

We analyze the model (^ by exploiting a finite-system 
density-matrix renormalization group (DMRG) method 
with open boundary conditions. ■^^'^^ Note that due to the 
two orbitals in one site, the number of bases for the sin- 
gle site is 16 and the size of the superblock Hilbert space 
grows as to^x16^, where m is the number of states kept 
for each block. Here, we treat each orbital as a site to 
reduce the size of the superblock Hilbert space to to^x4^. 
Then, the original zigzag chain is considered as a four-leg 
ladder and the finite-system sweep is performed along a 
one-dimensional path, as shown in Fig. 2. We note that 
in the superblock configuration, system and environment 
blocks are connected due to the on-site Coulomb inter- 
action as well as the inter-site electron hopping. Thus, 
we can increase m and accelerate the calculations due to 
the reduction of matrix size. In the present calculations, 
we keep up to m=300 states and the truncation error is 
kept around 10~^ or smaller. 

In order to determine uniquely the orbital state, it is 
useful to introduce an angle 9i to characterize the orbital 
shape at each site. We define new operators for electrons 
in the rotated frame as22i22i2i*S 



dipcr 



RiOO ( ^- 



(4) 






V4, Ch = tba = -V3t/4, ttb = 3i/4 (2) 



where the transformation matrix R{0i) is given by 



TABLE I: Angle 9i and the corresponding orbitals. 
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For reference, the correspondence of the angle Oi to the 
orbitals is summarized in Table I. 

Now the problem is how to determine the set of angles 
{6'i}. Naively thinking, a straightforward way is to find 
the set of {6{\ which minimizes the energy of the target 
state. However, due to the rotational symmetry in the 
orbital space, the energy does not depend on {6'i} in the 
electronic model including only Coulomb interactions i^ 
In such a situation, in order to determine {^i}, it is nec- 
essary to maximize the orbital correlations, as shown in 
Ref.f 

In general, to grasp what types of correlations develop, 
we investigate appropriate correlation functions. As for 
the orbital state, in the pseudo-spin representation for 
two orbitals, it is necessary to measure the orbital struc- 
ture factor, defined by 
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FIG. 3: DMRG results for orbital states and spin correlations 
in the ground state. The site dependence of {9i} representing 
orbital arrangement for (a) A=0 and (b) A=l. Spin correla- 
tion function measured from the center site of the lower chain 
for (c) A=0 and (d) A=l. (e) Electron densities in Zz^—r^, 
x^—y^, and optimal orbitals as a function of A. (f) Angle 6i 
at the center site of the system as a function of A. 



with f{'^Y.a0iaad-ia<T-dl0„di0a)/2, whcrc (• • •) denotes 
the expectation value. For each q, we find the set of {6'i} 
to maximize T(q). Among them, the optimal set of {^i} 
should be given by the one with the largest value of T(q). 
Namely, the optimal orbital at each site is determined so 
as to maximize the amount of the orbital correlation. We 
note that the evaluation of T(q) for any set of {6'i} can be 
carried out with just a single DMRG calculation by keep- 
ing appropriately the components of orbital correlations 
in terms of the original orbital basis set, due to the rela- 
tion f;- = T,- cos Oi + T^ sin 0, with T-^ = EA4aAaa - 



dlAb<r)/2 and T.^ 



EaidJaadiba 



- d^ rl- 



)/2. 



Although the set of {0i} cannot be determined from the 
energetical discussion of the two-orbital model, it is useful 
to consider the energy of the effective one-orbital model 
composed of the lower-energy orbitals specified by {6'i}. 
It has been found that the energy of the effective one- 
orbital model is actually minimized when the set of {^i} 
is determined so as to maximize the orbital correlation. 
Thus, we believe that our scheme to maximize the orbital 
correlation works well to find the optimal orbital state. 



III. GROUND STATE 

A. Spin and orbital structures 

Let us first clarify the orbital structure in the spin- 
singlet ground state, an important issue to understand a 
key role of orbital anisotropy in geometrically frustrated 
systems. Irrespective of A, we observe that T(q) be- 
comes maximum at q=0, indicating a FO structure. In 
Figs. 3(a) and 3(b), we show the optimal set of {6i} for 
A=0 and 1, respectively. The corresponding orbital ar- 
rangement in the bulk is also depicted. At A=0, two or- 
bitals are degenerate, but in general, orbital degeneracy 
is lifted to lower the energy by electron hopping, leading 
to an optimal path for the electron motion. Indeed, there 
appears a FO structure characterized by 6'i /7r~ 1.32, in- 
dicating that 3x^—r^ orbital is selectively occupied. The 
shape of the occupied orbital extends just along the x 
direction, not along the zigzag path, leading to a double- 
chain structure for the electron motion. Then, the spin 



exchange interaction along the zigzag path is expected to 
be suppressed in comparison with that along the double- 
chain direction. To clarify this point, we investigate the 
spin correlation function 



C,pi„(iJ) - {S?S? 



(7) 



with S'f =^^(/9irt— /5iTi)/2. In Fig. 3(c), the correlations 
measured from the center site of the lower chain is shown 
for A=0. We find that there exists an AFM correlation 
between intrachain sites in each of the double chain, while 
the spin correlation between the two chains is much weak. 
Thus, the zigzag chain is decoupled to a double-chain 
system due to the orbital anisotropy in terms of the spin 
correlation. 

On the other hand, the orbital anisotropy almost dis- 
appears at A=l, since electrons favorably occupy the 
lower 3z^—r^ orbital, which is isotropic in the xy plane, 
as shown in Fig. 3(b). In accordance with the variation of 
the orbital shape, electrons turn to move along the zigzag 
path as well as the double-chain direction. Then, the spin 
correlation appears between the two chains, as shown in 
Fig. 3(d). Note, however, that orbital degree of freedom 
is not perfectly quenched due to the level splitting. In 
fact, it is found that optimal orbitals at the edges deviate 
from those in the bulk. Namely, the orbital shape in the 
bulk extends to the double-chain direction in some de- 
gree, while it changes to become rather isotropic around 
the edges so that electrons can smoothly move near the 
edges. Thus, orbital degree of freedom is still effective, 
and there occurs the adjustment of the optimal orbital 
to the environmental lattice inhomogeneity. 

In order to gain insight into the change in the orbital 
state due to the level sphtting, we investigate the electron 
density in each orbital 



J2(P^r)/N, 



(8) 



as shown in Fig. 3(e). As A is increased, electrons are 
forced to be accommodated in the lower 3z'^—r^ orbital, 
but the electron density in each of 3z^—r^ and x^—ip or- 
bitals is found to change gradually without any anoma- 
lous behavior. Accordingly, the optimal a orbital state 
is smoothly transformed from 3x^— r^ to 3z^— r^ with in- 
creasing A, as shown in Fig. 3(f). On the other hand, we 
observe that optimal a orbitals are occupied, while (3 or- 
bitals are almost vacant, irrespective of A. This behavior 
is well understood by regarding the level splitting as the 
fictitious magnetic field in the orbital space. Namely, the 
direction of fully polarized orbital moment is rotated by 
the magnetic field. Thus, the present system is always 
regarded as a one-orbital system composed of optimal a 
orbitals, although we have considered the multi-orbital 
system. 



B. Effective spin model 

Let us here consider the effect of the appearance of an 
orbital arrangement on magnetic properties. For finite A, 
the lower Sz^— r^ orbital is occupied in the atomic limit, 
but strong electron-electron correlations cause a mixed 
orbital state of iz^—r^ and x^—ip orbitals. In fact, we 
have found that one optimal orbital becomes relevant due 
to the FO arrangement, while the relevant orbital is con- 
trolled by A. In order to understand magnetic properties 
of such a one-orbital system, it is useful to consider an 
effective model in the strong-coupling limit. For this pur- 
pose, we first express the Hamiltonian in terms of new 
operators by substituting Eq. Q as 



H ^ Y. t^^dl^A+^r'. + uY: 



PirfPiTi 



-U'Y,Pio.Pii3 -Aj2icoseift+smeifn,{9) 



where pi^^^dl^^d^^^, PiT=Y.a Pira, f-^ = T.M\aA^I3a + 

dla^diacr)/2, and the summation for the orbital index is 
taken over optimal a and /3 orbitals in the rotated frame. 
The hopping amplitude is given by 



tl 



fji^U 



■R{e,) 



W^'^y^uV^ 



-'i+ayji/r' 



(10) 



We note that U and U' (U^U') in Eq. @ are given by 
the same ones in the original Hamiltonian Eq. J^l, since 
the orbital basis set is transformed within the rotational 
invariant orbital space, and thus, the Coulomb integrals 
are expressed in the same form by Racah parameters^ 
Note also that, as mentioned in the previous subsection, 
the level splitting is regarded as the fictitious magnetic 
field for the orbital moment T? ^ while there appears the 
transverse field for T-^ . 

Magnetic and orbital ordering phenomena in multi- 
orbital systems have been actively discussed on the basis 
of an effective spin-orbital model, the so-called Kugel- 
Khomskii model^S*^ Here, we intend to obtain an under- 
standing of magnetic properties on the orbital-arranged 
background composed of the relevant a orbital. Thus, we 
consider a simplified spin model, by ignoring the irrele- 
vant j3 orbital. Then, we keep only the a-orbital parts in 
Eq. IjHJl, such as 

i,a,a" i 

-(A/2)^cos0iPi„. (11) 

i 

Now we apply the second-order perturbation with respect 
to electron hopping, and reach an effective Heisenberg 
spin model 

-ffAFM = Y. ^i.i+a(Si • Si+a - 1/4) + E,, (12) 
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FIG. 4: (a) Spin structure factor for A=0 and 1 in the ground 
state, (b) The A dependence of the peak position, estimated 
from the DMRG resuhs with Af=64. 



where Si is a spin-1/2 operator at site i, the AFM ex- 
change is given by 



J\M 



4|t 



l,i+a|2 



IV, 



(13) 



and £'s denotes the energy shift due to the level splitting, 
given by 



E, 



= -(A/2) ^ cos 6*1. 



(14) 



It should be noted that this effective model represents 
a ih%Tja% spin chain including nearest-neighbor exchange 
J\ along the zigzag path and next-nearest-neighbor ex- 
change J2 along the double chain. To discuss magnetic 
properties of the present eg-orbital model based on the 
effective model, let us here summarize the property of 
the zigzag spin chain, which has been revealed by inten- 
sive numerical and analytical investigationsAiSiiiiiSii^ In 
the zigzag spin chain, the ratio of competing exchanges 
J2/ Ji is a controlling parameter for the degree of spin 
frustration. At J2/J1—O, the system is equivalent to 
a single chain, which has a gapless spin-liquid ground 
state. With the increase of J2/J1, the spin frustra- 
tion becomes effective, and the ground state is changed 
from a critical spin-liquid phase to a gapped dimer phase 
at J2/Ji=0.24>i^ The dimer phase continues up to the 
limit of J2/ Ji=oo, while the spin gap decreases expo- 
nentially,^^ since the system turns to be independent two 
chains for J2/ Ji—00. On the other hand, the spin cor- 
relation exhibits a commensurate AFM peak at q—n for 
0< J2/ Ji<l/2, while the peak changes into an incommen- 
surate one for J2/Ji>l/2, and the peak position gradu- 
ally moves to q—TT/2 for infinite J2/ Ji^^^ 

In order to confirm that magnetic properties are actu- 
ally described by the zigzag spin chain, let us discuss the 
behavior of the spin correlation. In Fig. 4(a), we show 
our DMRG results of the spin structure factor 
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(15) 



with <5f =X)r('°i^T^PiTi)/2- Note that the results are de- 
picted with a scalar wave number, which is defined by re- 
garding the zigzag path as a linear chain, for comparison 



with the characteristics of the zigzag spin chain described 
above. For A=0, we clearly observe a peak at g=7r/2, 
consistent with that of the zigzag spin chain with large 
J2I Jii since J2/Ji=64^ for the "ix^—r^ FO arrangement. 
Note that we evaluate J\ and J2 at the center of the 
chain, by ignoring the site dependence in Ji.i+a- The ef- 
fect of the site dependence of Ji,i+a on the spin gap will 
be discussed in the next section. In Fig. 4(b), we show 
the A dependence of the peak position of the spin struc- 
ture factor. With increasing A, the peak position changes 
away from g=7r/2 towards q=vr, as expected by analogy 
with the zigzag spin chain, since J2/ Ji approaches unity 
such as J2/Ji = 1.61 for A=l, J2/Ji = 1.24 for A=2, and 
>/2/>/i=l-08 for A=5. Note that in the zigzag spin chain, 
the peak locates at g~0.56l7r for J2/ Ji=l.~ In fact, the 
peak in the present Cg-orbital model approaches this po- 
sition with increasing A, as shown in Fig. 4(b). Thus, the 
behavior of the spin structure factor is well understood 
by analogy with the zigzag spin chain. 



IV. SPIN EXCITED STATE 

Now we turn our attention to the spin excitation. In 
order to discuss the spin gap, in general, we consider the 
lowest-energy state among spin-triplet states. One may 
expect that such a triplet state is also described by the 
Heisenberg model, but we should note that the change of 
the spin state influences on the orbital state as well. Such 
effect is already found in a two-site problem, as shown in 
Fig. 5(a). In the upper panel, we depict an AFM/FO 
configuration in the spin-singlet ground state. On the 
other hand, in the spin-triplet excited state, an antiferro- 
orbital configuration turns to be favorable to gain kinetic 
energy by electron hopping due to the Pauli's exclusion 
principle, as shown in the lower panel. 

In order to clarify the characteristics of the orbital ro- 
tation due to the spin-triplet excitation, we investigate 
the spin-orbital correlation function, defined by 



Cso(k,l;i,j) = ((3/4 + Sk 



Si)i;^r/ 



(16) 



which represents the orbital correlation on the condition 
that the spin state of bond (k, 1) is triplet. In Fig. 5(b), 
we show the spin-orbital correlation function in the spin- 
triplet excited state. Here, we set i=k to consider the 
orbital rotation around the spin-triplet bond. We ob- 
serve that staggered orbital correlation is enhanced on 
the spin-triplet bond, as was expected from the two-site 
problem, while the staggered component rapidly decays 
for long distances from the spin-triplet bond. Namely, 
the orbital rotation occurs only close to the spin-triplet 
bond. Note that orbital correlations near the edges are 
found to increase due to the deformation of the optimal 
orbital around the edges, as shown below. 

Since the spin-triplet pair moves in the system with 
carrying the orbital rotation, it is expected that the 
averaged orbital structure should be different between 
the spin-singlet ground state and the spin-triplet excited 
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FIG. 5: (a) Schematic view of spin-orbital configuration, (b) 
Spin-orbital correlation function in the spin-excited state, (c) 
The optimal set of {9i} in the ground and spin-excited states, 
(d) The system-size dependence of the difference in the energy 
shift between ground and spin-excited states. 
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state. In fact, in the excited state, a FO arrangement is 
also found to appear and the system is regarded as a one- 
orbital system in a similar way to the case of the ground 
state. However, as shown in Fig. 5(c), we clearly observe 
that optimal orbitals in the excited state are changed 
from those in the ground state. This behavior is under- 
stood from a viewpoint of the adjustment of the optimal 
orbital. Namely, the orbital state is reorganized so as 
to provide the most appropriate orbital arrangement for 
the spin state. Note that in the excited state, the orbital 
shape around the edges becomes isotropic in comparison 
with those in the bulk similarly to the case of the ground 
state. 

In order to ascertain this characteristic change of the 
orbital arrangement, we evaluate the difference in the 
energy shift. 



SK = Ei'^-Ei'K 



(17) 



where Eg and Eg denote the values of Eg in the ground 
and excited states, respectively. In Fig. 5(d), we show the 
system-size dependence of SEg. It is observed that SEg 
grows and converges to a finite value with increasing the 
system size, suggesting that SEg will contribute to the 
spin gap. Concerning exchange interactions { Ji,i+a}, we 
obtain the ratio at the center of the chain as J2/ Ji=1.61 
in the ground state, while J2/ Ji = 1.62 in the excited state 
for A=l from the DMRG results of A^=40. One may con- 
sider that the difference seems small, but as we will see 
later, even such small difference leads to the significant 
change in the spin gap. In short, the present system can 
be regarded as a one-orbital system even in the excited 
state, but the variation of the orbital arrangement ap- 



pears in the changes of Eg and { Ji,i+a} in the effective 
Heisenberg model. Thus, spin and orbital degrees of free- 
dom dynamically interact with each other, and the spin 
excitation involves the changes of not only spin state it- 
self, but also orbital state. 

Now let us discuss the behavior of the spin gap at A=l. 
From the estimation of J2/ Ji=1.61 in the ground state, 
the ground state is expected to correspond to the gapped 
dimer phase of the zigzag spin chain iii*iSiii To clarify this 
point, we investigate the system-size dependence of the 
spin gap 

A,pi„ = EiN/2 + 1, N/2 - 1) - E{N/2, N/2) (18) 

up to 7V=64, where E{N^,Ni) denotes the lowest energy 
in the subspace with A'^^ up- and N^ down-spin electrons. 
As shown in Fig. 6(a), the spin gap seems to converge to a 
finite value in the thermodynamic limit. For comparison 
with the zigzag spin chain, we also evaluate the spin gap 
for two types of AFM Heisenberg models. First, we take 
a uniform model, assuming that AFM exchanges are uni- 
formly given by Ji and J2, which are estimated at the 
center of the chain for the ground state. In Fig. 6(a), 
such results are denoted by crosses. Corresponding to 
the gapped dimer phase, we observe a finite spin gap in 
the thermodynamic limit, but significant deviation exists 
between the results of the eg-orbital Hubbard model and 
those of the uniform Heisenberg model. 

It is emphasized here that such deviation indicates the 
effect of the orbital-arranged background which should 
be different between ground and excited states. Then, 
we attempt to re-estimate the spin gap by considering 



the following two points: (i) The optimal orbital is not 
uniform in the system, especially around the edges, as 
shown in Fig. 3(b). (ii) The orbital arrangement is found 
to be different between ground and excited states. To 
take account of these points, we deal faithfully with the 
effective Heisenberg model Eq. H12|) with the optimal set 
of {6i} for both ground and excited states. The variation 
of { Ji4+a} causes the change in the energy of the spin 
system. There also occurs the difference in the energy 
shift due to the level splitting, as shown in Fig. 5(d). 
Note that we calculate the spin gap for the model with 
site-dependent exchange interactions by exploiting the 
DMRG algorithm for random spin systems^ The results 
are denoted by open circles in Fig. 6(a). We find that 
the spin gap is appropriately described by the effective 
Heisenberg model with the optimal orbital rather than 
the simple uniform Heisenberg model. 

Although some deviation still remains, it is improved 
by increasing m to keep more states in the renormal- 
ization step for the Cg-orbital Hubbard model. Another 
aspect of the remaining deviation would be the effect of 
charge and/or orbital fluctuations due to the smallness of 
the Coulomb interaction to arrive at the strong-coupling 
limit. We note that larger energy is obtained for the effec- 
tive Heisenberg model than that for the Cg-orbital Hub- 
bard model in both ground and excited states, since the 
fluctuation effects are ignored. However, such increase of 
the energy is found to be different between ground and 
excited states. As a result, the spin gap of the effective 
Heisenberg model becomes small in comparison with the 
e^-orbital Hubbard model. To see the fluctuation effects 
on the spin gap, we show the A and U' dependence of 
the spin gap in Fig. 6(b) and 6(c), respectively. With 
increasing A, J2/J1 approaches unity and the spin gap 
increases, as expected by analogy with the zigzag spin 
chaiufii while the deviation remains significant due to 
the effect of charge fluctuation. On the other hand, as 
shown in Fig. 6(c), the deviation is reduced with increas- 
ing U' , indicating that the spin gap is well reproduced by 
the effective Heisenberg model with the optimal orbital, 
when the charge and orbital fluctuations are suppressed. 
Thus, we conclude that active orbital degree of freedom 
influences on the spin gap, suggesting that the spin gap 
should be considered as a spin-orbital gap. 



V. SUMMARY 

In this paper, we have discussed ground-state proper- 
ties and spin gap of the Cg-orbital Hubbard model on the 



zigzag chain. We have found that the system is regarded 
as an effective spin system on the orbital-arranged back- 
ground due to the orbital arrangement. At A=0, two 
orbitals are degenerate, but 3x^—r^ orbital is selectively 
occupied so as to suppress the spin frustration due to 
the orbital anisotropy. On the other hand, the orbital 
anisotropy is controlled by A. With increasing A, the 
lower 3z^—r^ orbital is favorably occupied and the or- 
bital anisotropy almost disappears, leading to the revival 
of the spin frustration. Note that the lower 3z^—r^ or- 
bital is also occupied in the atomic limit C/'— >oo for A>0. 
Then, the system is corresponding to the gapped dimer 
phase of the zigzag spin chain. In actual compounds, A 
and U' could be controlled by hydrostatic or chemical 
pressure, indicating the control of the degree of the spin 
frustration and a possible realization of the gapped phase 
of the zigzag spin chain. 

However, orbital degree of freedom is not completely 
quenched due to the level splitting even at A=l. In fact, 
there appears a spatially non-uniform pattern of the op- 
timal orbital especially around the edges so as to adjust 
to the lattice inhomogeneity In addition, the orbital ar- 
rangement is flexibly deformed according to the change 
of the spin state. In order to discuss the spin-gap for- 
mation in multi-orbital systems, it is quite important to 
consider the interplay of spin and orbital degrees of free- 
dom, even when orbital degree of freedom is suppressed 
due to the level splitting. The spin gap is essentially de- 
termined by the spin-orbital coupled excitation. Thus, it 
is quite natural to consider the spin gap as a spin-orbital 
gap. 
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